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ABSTRACT 



The interrelation of hard X-ray (HXR) emitting sources and the underlying physics of electron acceleration and transport presents one 
of the major questions in high-energy solar flare physics. Spatially resolved observations of solar flares often demonstrate the pres- 
ence of well-separated sources of bremsstrahlung emission, so-called coronal and footpoint sources. Using spatially resolved X-ray 
observations by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and recently improved imaging techniques, 
we investigate in detail the spatially resolved electron distributions in a few well-observed solar flares. The selected flares can be 
interpreted as having a standard geometry with chromospheric HXR footpoint sources related to thick-target X-ray emission and the 
coronal sources characterised by a combination of thermal and thin-target bremsstrahlung. Using imaging spectroscopy techniques, 
we deduce the characteristic electron rates and spectral indices required to explain the coronal and footpoint X-ray sources. We found 
that, during the impulsive phase, the electron rate at the looptop is several times (a factor of 1 .7 - 8) higher than at the footpoints. The 
results suggest that a sufficient number of electrons accelerated in the looptop explain the precipitation into the footpoints and imply 
that electrons accumulate in the looptop. We discuss these results in terms of magnetic trapping, pitch-angle scattering, and injection 
properties. Our conclusion is that the accelerated electrons must be subject to magnetic trapping and/or pitch-angle scattering, keep- 
ing a fraction of the population trapped inside the coronal loops. These findings put strong constraints on the particle transport in the 
coronal source and provide quantitative limits on deka-keV electron trapping/scattering in the coronal source. 
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1. Introduction 

The X-ray images of solar flares with spectral capabilities 
taken by the Hard X -ray Telescope (HXT) on board of Yohkoh 
dKosugi et alJ 1 1 99 lh and Reuve n Ramaty High Energy Solar 
Spectroscopic Imager (RHESSI) (iLin et al.ll2002l) have revealed 
a wealth of information on spatially resolved distributions 
of energetic electrons in solar flares (see e.g . iDennis et al.l 
1201 It iHolman et all 1201 it iKontaret al.l 1201 laL as recent re- 
views.) These spatially resolved observations indicate ap- 
pearances of distinct coronal and footpoint X-ray sources. 
Hard X-ray (HXR) footpoint sources are usually descr i bed in 



~ 10-30 keV and are usually due to thermal and non- 
thermal bremsstrahlung, which indicates plasma he ating and the 
prese nce of non-thermal electron populations (e.g. lEmslie et al.l 
120031) . Spectroscopically, there is no clear boundary between 
the two processes, with the spectrum transiting smoothly from 
thin-target emission in the corona to thick-target emission within 
the range 10-30 keV. The looptop sources are believed to 
be connected to the process of energy re lease in flares (e. g. 
iPriest & Forbesll2002l: IHolman et al.ll201 lh . Some models place 
the magnetic reconnection at the looptop, where the X-ray 
emission is a consequence of the plasma heating and particle 
acceleration ( Jakimiec et al. 19981: Petrpsian & Donaghvlll999l: 



terms of collisional thick- t arget models (iDubov 19631 



1969t lArnoldv et all [19681 iLingenfelterl Il969l IBrownl 1 197 it 
Svroyatskii & Shmeleval Il972t iBrownl Il973al Lin & Hudson 



19761) . where the electrons accelerated in the corona travel 



downwards inside magnetic loops marking footpoint sources 
in dense chromospheric pl asma by emission in X-rays, EUV, 
and white-light emissions (iPhillips et al.ll2008h . The observa- 
tions show that the HXR source heigh t in the chromosphere 
decreases with increasing X-ray energy ([Matsushita et al.lll992t 
Aschwanden et al.l 120021: iKontar et all l2Q08h ISaint-Hilaire et all 
20 lOt iBattaglia & Kontd 1201 1 ah as anticipated from down- 
ward propagating electron beams. In addition, the characteristic 
sizes of HXR sources decrease with increasing energy, which 
is consistent with collisi onal electron transport along converg- 
ing magnetic flux tubes (IKontaret al]l2008t IBattaglia & Kontan 
l20lTllFedun^ et al]l2012tlXu et alJl2012l) . 

Coronal X-ray sources are prominently visible at and below 



iBastian et alj|2007tlKontar et alj|201 lbtlGuo et alj|2012h . Other 
models propose that the X-point is placed higher in the corona, 
and the looptop X-ray source shows that the interaction between 
a do wnward reconnection flow and plasma fills a soft X-ray 
loop dTsuneta et al.lll997h . The recent observations of energy- 
dependent HXR source sizes inside the dense coronal loops 
(IXu et al.ll20Tj8tlKontar et al.ll201 lbtlBian et alJl20Tll:lGuo et ail 
l2012r.lTorre et al.l2012h also suggest within-the-loop electron ac- 
celeration and transport. 

iPetrosian et al.l d2002l) presented a comprehensive analysis 
of 18 Yohkoh events, finding that the spectral index of the 
looptop source is softer than the footpoints on the average 
by about 1, although this result is limited by the poor energy 
resolution of the Yo hkoh images, with only four energy ranges 
dKosugi et al.ll99lh . With the launch of RHESSI, detailed imag- 
ing spectroscopy with e nergy resolut i on ~ 1 ke V has become 
widely available (e.g. lEmslie et al.l 120031 : IBattaglia & Benzl 
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1200ft iPiana etafll2007t iKundu et al]l200l) . Thus, Emslie et al. 
(120031) found that for the GOES X-class 2002 July 23 flare the 
coronal source is consistent with the thermal source of around 
4 x 10 7 K and the footpoint-like sou rces that have power-law 
spectra with different spectral indices. Bat taglia & Bend d2006l) 
presented a systematic study of the relation between looptop 
and footpoint sources for a few flares observed by RHESSI. 
They found soft-hard-soft behaviour in both the coronal source 
and the footpoints. Assuming the scenario where the same 
electron population produces the coronal source by thin-target 
bremsstrahlung and the chromospheric source by thick-target 
bremsstrahlung, one would expect a difference of ^ 2 between 
the two photon spectral indices. They found that the average 
of all mean differences is =* 1.8. However, the difference is 
considerably larger than 2 for two out of five events, and 
smaller than 2 for the other three events analyzed. The coro- 
nal source is nearly always softer than the footpoints. The 
footpoint spectra differ significantly only in one event out of 
five. While the observations are consistent with acceleration 
in the coronal source and a subsequent propagation into the 
chromosphere, their result excludes the simplest scenario of 
electrons free-streaming from the acceleration region towards 
the footpoints, which indicates that other transport processes 
must be considered, e. g. wave-partic le interaction can accoun t 
for spectral differenc e s large r than 2 (lHannah & KontaJl201 ll) . 
iTomczak & Ciborskil d2007l) used HXT/Yohkoh data of 117 
flares to investigate the ratio of photon counts from looptop 
to footpoint sources, finding a correlation between this ratio 
and height and attributing the effect of loop convergence in 
80% of the flare sample to inferred mirror ratios lower than 
2.1. Coronal trapping is a natural consequence of the expected 
convergence of the magnetic field between the corona and the 
chromosphere. Trap- plus-precipitation models have been pro- 
posed and developed dTakakura & Kail 1966^ Melrose & Brownl 



19831: iMacKinnonl 1198 



197ft iLea ch & Petrosi anl 11981 

Alexander! 1 990t iMcClementsll 1 9901; IMacKinnonl 1 99 It [Fletche 
1997t iFletcher & Martensl Il998l Fark & Fleishmanl 1201 Oh 



Also, considerabl e observation a l evid e nce of trapping has 

been provided ( Ramaty et al.l 11994 ; Alexander & Metcalfl 

1997t lAschwanden et all 119981: iMetcalf & Alexander! | 1999[ 
Lee et al. " p000t lAschwanden et al.l 119991: iLee & Garvl bood 
Leeetal.l 120021: I Alexander & Metcalfl 120021 IMelnikov etall 



2002t iFleishman et al.l 1201 lb iHuang & Lil [201 ll). Analysing 
microwave brightness maps, IMelnikov et al. ( 20021) applied a 
collisionless adiabatic trap model to explain looptop sources 
at 17 and 34 GHz observed in four flares. They concluded 
that the electron distribution along the loops must be highly 
inhomogeneous, with the energetic electrons concentrated at the 
upper part of the loop. On the other hand, microwave loop-like 
brightness structures of different flares could only be reconciled 
with models of uniform magnetic field and h omogeneous 
spatia l distribution of no n-thermal electrons (IKundu et alJ 
l200l . iNindos et al.l d2000l) successfully reproduced the main 
spatial morphology of 5 and 15 GHz brightness maps of a weak 
flare using a non-uniform magnetic field model and uniform 
distribution of electrons, in agreement with model simulations 
dAlissandrakis & Preka-Papademal[T984l iKlein & Trotted fl98l 
ISimSes & Costal 20061 l2010l) . Those examples evidence how 
trapping conditions can change from flare to flare, influencing 
the electron distribution along flaring loops. Moreover, it 
seems that the microwave emission maps indirectly show the 
trapped population, as the emission also depends strongly on 
the magnetic strength. However, HXR images will mainly 
show the precipitated population as footpoint sources, only 



eventually showing a non-thermal HXR coronal source. As the 
collisions are energy-dependent, one can expect that deka-kev 
HXR-producing electrons and MeV microwave-producing 
electrons will evolve in slightly different ways during a solar 
flare (e.g. lLee et al 1 l2002tlGimenez" de Castro et al 112012b . 

Importantly, electron scattering (often str ong) is required to 
achieve efficient acceleration of particles (e.g. iMiller et"ai1ll997l 
iMillerlll 9981: iPetrosianll 9991: iBian et al.ll2012l) and to explain the 
lack of strong HXR anisotropy in observations, while the trans- 
port of particles affects how the particles are accelerated (see e.g. 
iBvkov & Fleishmanll2009l:lBian et al J|20 1 2t lPetrosianll20 1 2h . 

Despite the importance of comparative studies between the 
coronal sources and footpoints, there were no systematic com- 
parisons of the number of electrons in the footpoints and the 
coronal sources for the events with visible footpoints. We em- 
phasize that the inference of imaging spectroscopy parame- 
ters relies strongly on the precise knowledge of the X-ray 
source area, which was only developed fairly recently (e.g . 



Schmahl et al 1120071: lKontaretalJl2008l: I Dennis & Pernakll2009l: 



Kontar etaU l2010t) . In addition, the full Spectral Response 



Matrix (SRM) of RHESSI detectors (including non-diagonal 
terms) for imaging spectroscopy became available relatively re- 
cently, since 2006 February. 

We present a comparative analysis of the electron distribu- 
tions in footpoints and coronal sources for four well-observed 
flares. We find for the first time the electron rate required to ex- 
plain the coronal and chromospheric X-ray emissions. In Sect. 
[2] we present the observations and in Sect. [3] the methodology 
to retrieve the electron rates from imaging spectroscopy. Our re- 
sults are presented in Sect.|4]and discussed in Sect. [5] The work 
is summarized in Sect. [6] 

2. Observational data 

We selected four well-observed events in which the common 
loop structure could be identified: two footpoint sources at 
higher photon energies and a looptop source at lower energies. 
We avoided events where looptops and footpoints overlap due to 
the projection of the loop geometry. Also, we only selected flares 
where a non-thermal source at the looptop could be resolved. 
The flare list with details is presented in Table [TJ Different 
aspects of the selected fla res were previously investigated by 
oth er authors: 2002 Ju ly 23 (lEmslie et al.l2003l). 2003 November 
02 dSilva et al.l 120071) . 2011 February 24 dBattaglia & Kontar] 
1201 lbh . 



3. Methodology 

3. 1 . Imaging spectroscopy 

For each of the fou r flares, we constructed CLEAN images 
dHurford et al.ll2002l) with a pixel size of 1", using front detec- 
tors 3 to 8, for 19 logarithmically binned energy bands from 10 
to 100 keV. Since we rely on the source sizes taken from the 
reconstructed images, we verified the b est CLEAN beam width 
for each flare. iDennis & Pernakl (120091) and iKontar et al l (1201 Ol) 
have pointed out that CLEAN images usually have systemati- 
cally larger sizes than other algorithms when using the default 
beam width factor of 1.0. To ensure the best possible determi- 
nation of the source sizes using CL EAN images, we app lied the 
visibility forward-fitting procedure (ISchmah l et al. 2007) on the 
footpoints of each flare, adjusted the CLEAN beam size, and 
re-calculated the images until the FWHM of the footpoint from 
both algorithms had the same size. The integration time interval 
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Table 1. Selected solar flares: date and time used for analysis, GOES peak and classification events, RHESSI spectral fitting results, 
source sizes measured from HXR images, and electron rate results for looptop Nlt and footpoint Nfp sources. 



Flare 




A 


B 


C 


D 


Date 




2002 July 23 


2003 Nov 02 


2011 Feb 24 


2011 Sept 24 


Time 


UT 


00:27:26 


17:15:54 


07:29:40 


09:35:53 


At 


sec 


284 


246 


176 


82 


GOES peak 




X5.1 


X8.3 


M3.5 


X1.9 


GOES time 


UT 


00:36 


17:24 


07:35 


09:40 


CLEAN beam 




2.3 


2.0 


1.9 


2.4 


Imaging spectroscopy fitting results 












< hVF > 


lO^cnrV 1 


77 ± 12 


169+12 


0.7 ±0.1 


7±2 


Slt 




3.7 ± 0.2 


2.9 ±0.1 


3.3 ± 0.2 


3.6 ±0.4 


Nfp 


10 3 V 


32 ±5 


153 ± 11 


1.1 ±0.1 


5.2 ± 1.1 


S F p 




4.0 ±0.1 


3.9 ±0.1 


3.8 ±0.1 


3.8 ±0.1 


Spatially integrated spectral fitting results 












EM 


10 49 cnT 3 


3.8 ± 0.4 


4.2 ±0.6 


0.24 ± 0.02 


0.8 ±0.1 


T 


MK 


33 


34 


21 


34 


N 


10 3 V 


27 ±3 


122 ± 11 


2.0 ±0.1 


8.9 ± 1.0 


d 




3.9 ± 0.1 


4.0 ±0.1 


4.0 ±0.1 


3.8 ±0.1 


h 


1 r\ 11 _7 

10"cm 1 


2.7 


1.5 


0.6 


1.9 


Characteristic sizes of the sources 












A 


, A |S 1 

10 cm" 


0.6 


2.2 


0.9 


0.3 


F\ 
U 


1U cm 




o.U 


O.J 


0.4 


V* 


10 27 cm 3 


0.5 


1.8 


0.8 


0.2 


L 


10 8 cm 


5.2 


4.3 


8.3 


4.7 


Electron rates 












Nlt 


10 3 V 


54+11 


258 ± 38 


1.6 ±0.3 


8±3 


N F p (fully ionized FP, without albedo) 


10 3 V» 


32 ±5 


153 ± 11 


1.1 ±0.1 


5.2 ± 1.1 


NltINfp (as above) 




1.7 + 0.6 


1.7 ±0.4 


1.5 ±0.4 


1.4 ±0.9 


N F p (neutral FP, with albedo) 


10 3 V 


6± 1 


27 ±2 


0.4 ±0.1 


0.8 ±0.2 


NltINfp (as above) 




9.3 ± 3.0 


10.1 ±2.3 


4.5 ± 1.4 


7.8 ±5.2 



during the impulsive phase and the CLEAN beam used for each 
event are presented in Table [TJ The CLEAN beam factor found 
for these flares lay in the range 1.9 - 2.4, suggesting the opti - 
mum value is ^ 2 (iDennis & Pernakll2009l: IKontar et al.l[2010h . 
instead of the default value of 1. Fig. [TJ shows CLEAN maps 
for the four flares at lowest (10-11.3 keV) and highest (88.6-100 
keV) energy bins. 



3.1 .1 . Regions of interest selection 

Within each image, we selected regions of interest (ROI) which 
capture looptop (LT) and footpoint (FP) sources for fu rther 
imaging spectroscopy using OSPEX dSchwartz et al.ll2002b . For 
the flares considered here, two f ootpoints can be easi l y identified 
and no ted i n the literature, e. g. lEmslie et al.l (l2003l) . ISilva et al.l 
d2007l) . and lBattaglia & Kontaj d201 lbl) for flares A, B, and C, 
respectively. The looptop ROI was chosen to include the main 
emitting regions at low energies, but without being too large, to 
avoid contribution from the FPs. The looptop and footpoint ROIs 
for each flare are shown in Fig. [TJ The background was estimated 
taking the sum of pixels outside the ROIs in each image, then 
normali zed for the area of each R OI, similar to an approach ap- 
plied bv lBattaglia & Benzl (I2006I) . This method probably overes- 
timates the background noise, especially at lower energies; nev- 
ertheless, it is useful to identify the reliable energy ranges for the 
fitting procedure. 



3.1.2. Spectral fitting 

While the detailed analysis of HXR spectrum in some events, 
e. g. 2002 July 23 flare, demonstrated that the actual spectrum 



does contain deviations fr om a simple power-law dPiana et al.l 
I2003I: IKontar et al.ll2011al) . we assume a single power-law in 
the electron spectrum for the purposes of our analysis. For each 
flare, the looptop source was fitted with a thin-target model and 
the footpoint source s were fitted using a thick-target model (see 
IKontar et al.ll201 lal for details). Assuming looptop source is a 
thin-target emission, the photon flux detected at the Earth is the 
convolution of the electron distribution and the bremsstrahlung 

cross-section. 

Following lBrownl d 1 97 ll) and lBrown et al.1 d2003l ). we define 
the mean target proton density h = V -1 J n(r) dV, where n(v) is 
the plasma density, V is the volume of the LT emitting region. 
We can then write the photon spectrum at distance the Earth as 



hrie) 



i r 

— ^ <n 

47TR 2 Je 



VF(E) > Q(e, E) dE, 



(1) 



where the bremsstrahlung cross-section Q(e , E) following dHaud 
119971) and F(E) is the meciYi electron flux spectrum defined by 



F(E) = — F(E,r)n(r)dV. 
nV J v 



(2) 



From the OSPEX dSchwartz et al.ll2002l) fit of the photon flux 
from looptop ROI with the model given by Eq. ([TJ, one finds 
< hVF{E) > without any additional assumptions, where 



< hVF(E) >=< hVF > 



, E>E 



(3) 



normalised to < nVFo > [electrons cm 2 s '], which is the prod- 
uct of the mean target density h and the energy-integrated (above 
the low energy cut-off Eq) mean electron flux Fq — J F(E)dE 
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Fig. 1. CLEAN maps for each flare: (a) 2002 July 23, (b) 2003 November 02, (c) 201 1 February 24, and (d) 201 1 September 24, 
showing the emission at 10-11.3 keV (black contours) and 88.6-100keV (blue contours) at 30, 50, 70, 90% of maximum of each 
energy bin. The ROIs for the imaging spectroscopy are defined in gray lines: footpoints (continuous lines) and looptop (dashed 
lines). Integration times are shown in each frame and in Table Q] The CLEAN beam width factor used to construct each map 
according to the best agreement with visibility forward-fitting is also presented in Table [T] For each, a total of 19 maps were made 
with energy bins logarithmically spaced between 10 and 100 keV. 



in the volume V. The energy-integrated electron rate Nu [elec- 
trons s~'] required to explain the observed thin target emission 
in the looptop is simply 



where S is the cross-section area of the loop. We can estimate 
F by 



Nlt = FoS, 



(4) 



< nVF > < nVFo > 



FiV 



nLS 



(5) 
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where L is the length of the looptop source in the direction along 
the loop (hence V — LS). Substituting Eq.[5]into Eq.|4]we have 



N 



LT 



< nVF > 
hL 



(6) 



The thermal plasma density h and the size L are obtained from 
observation data (see Sect. 13.1. 31 . 

For the footpoint emission, we assume a thick-target model, 
so that the observed photon flux at the distance R is given: 

i r°° e l e \ Sfp ~ 2 

^^I^yUJ Q ^ E)dE > (7) 

where Nfp is the electron precipitation rate at footpoints [elec- 
trons s _1 ], and K = 2ne A K, e is the electronic charge and A is 
the Coulomb logarithm (assumed =* 22.7 for fully ionized chro- 
mosphere and 7.1 for neutral chromosphere, see Sect. 14.1) . Here 
we emphasize that 5fp and Nfp are the spectral index and the 
electron precipitation rate of the electrons entering into the foot- 
points. If the electrons are unaffected by the transport within the 
loop, this spectral index should match 5lt an d the electron rate 
that of looptop source Nlt- Both fi ts included a ther mal compo- 
nent, using full CHIANTI model dDere et alJl2009h . with con- 
tinuum and spectral lines. A spatially integrated spectrum was 
also made for each flare using the same time intervals and fit- 
ted with a thermal plus thick-target component. The thick-target 
model was set as a single power-law. The fitting results (Eqs. 
Q] and Q for the imaging and spatially integrated spectroscopy 
are presented in Table Q] The photon spectra of the LT, FP, and 
spatially integrated sources are plotted in Figs. [2] to E] for each 
flare, respectively, showing the fitting components (thermal and 
non-thermal) and fitting residuals. 



3.1.3. Image-derived parameters 

The overall geometry of the selected flares can be interpreted as 
a single loop-like magnetic structure with two footpoints and a 
coronal source at the apex of the loop. We defined a loop base- 
line as the line connecting the two footpoints, which is used as a 
reference to measure the sizes of the sources at the looptop. The 
thermal plasma density can be calculated using the observations 
of thermal X-ray emission. The isothermal fit to the spatially in- 
tegrated RHESSI spectrum (right-hand frames in Figs. [2] to [5} 
gives emission measure EM — n 2 V,i„ where the volume of the 
thermal source can be estimated as V,h = AD, A is the area of the 
source, D is the cross-section diameter of the source, and assum- 
ing uniform plasma density h within the volume V,h . Using the 
emission maps at 10-11.3 keV, D was estimated by measuring 
the length of this source at 50% level in the direction orthogonal 
to the baseline and A as the source area at 50% level of the max- 
imum emission. Using those measurements, the plasma density 
is then obtained by n = {EM I DA) 1 1 2 and directly inserted into 
Eq.|6l resulting in 



Nlt = 



/ DA \ 1/2 < hVF 



\EMI 



(8) 



Similar results for EM are obtained using the imaging spec- IKontar etal 
troscopy and selecting looptop regions only (left-hand frames 
of Figs. [2]to|5]l, but the imaging spectroscopy gives larger uncer- 
tainties. The values found for the plasma density «/10 n cm -3 
are 2.7, 1.5, 0.6, and 1.9 for flares A, B, C, and D, respectively. 
The geometry aspects of the flares and plasma density values 
inferred are presented in Table [I] 



We define the length L of the non-thermal source at the loop- 
top as the size (at 50% of maximum emission) in the direction 
parallel to the footpoint baseline and crossing the point with 
the maximum emission. L is measured at a sufficiently high en- 
ergy, where the contribution from thermal emission is not sig- 
nificant and where the source is spatially resolved. The non- 
thermal looptop sources can be easily identified at 42.8 - 48.3 
keV, 42.8 - 48.3 keV, 20.7 - 23.4 keV, and 29.8 - 33.6 keV for 
each flare, respectively (Pig. |6j»; the measured lengths L for each 
flare are presented in Table [1] 



4. Observational results and analysis 

Using Eq. [8] with the values found above, the electron rate at 
the looptop Nlt can be calculated and is summarized in Table 
\T\ along with the electron rate at the footpoints obtained using 
thick-target forward fit for the photon flux spectrum. 

Spectral index analysis of these four flares reveals that all 
flares have a similar footpoint spectral index, which is close to 
3.9 (see Table[TJ. The differences between the spectral index are 
6 LT - S F p = -0.3 + 0.2, -1.0 + 0.1, -0.5 + 0.2, and -0.2 + 0.4 
for each flare. This app ears consistent with previous findings, 
e.g. lEmslie et al.l (120031) reported that the photon spectral in- 
dex difference is less t han 2 for the 2002 July 23 flare and 
iBattaglia & Benzl (|2006j) reports a fairly broad range of spectral 
index differences between footpoints and coronal sources. 

In general, our results suggest that in all flares the precip- 
itation rate into footpoints is smaller than the electron rate re- 
quired to explain the coronal HXR emission, i.e. Nlt /Nfp > 
1. Assuming fully ionized chromospheric footpoints and no 
albedo, Nlt INfp for the four flares appears around 1.6, which is 
slightly larger than unity. This is expected in an idealized model, 
where non-thermal electrons propagate without interaction, so 
that Nlt — Nfp- However, as discussed in the next section, the 
effects of neutral chromosphere and albedo can increase this ra- 
tio by a factor of ~ 2 - 6. 

4. 1 . Effects of photospheric albedo and neutral 
chromosphere 

The solar atmosphere above the regions where the X-ray 
bremsstrahlung is produced during flares is optically thin, while 
the dense photospheric layers below those regions are opti- 
cally thick to X-ray photons. As a consequence, photons emitted 
downwards can be efficiently Compton backscattered by atomic 
electrons in the photosphere dTomblinl 19721 Santangelo et al.l 
119731: iBai & Ramatvlll978l: iMagdziarz & Zdziarskilll995l) . The 
observed X-ray spectrum from flares is then a combination of 
the primary and albedo photons. The reflectivity is dependent 
both on the primary photon spectrum and the flare position on 
the Sun relative to the observer, or more precisely, the helio- 
centric angle of the source. The backscattered radiation can be 
significant and modify the primary X-ray spectrum. More im- 
portantly for our analysis here, when the albedo contribution is 
taken into account, the number of e lectrons required to explain 
the observed phot on flux changes (lAlexander & Brownl 120021: 



2006). New fittings taking the albedo correction 



(IKontar et al.l l2006) into account give expectedly lower numbers 
for the electron rate at footpoints Nfp/10 35 electrons s~ l : 18 + 3, 
69 + 5, 0.9 + 0.1, 2.1 + 0.5, or factors of 1.8, 2.2, 1.2, 2.5 lower 



than the electron rate found without considering the albedo pho- 
tons. The effect is more pronounced for flares closer to the centre 
of the solar disk, while less significant for limb flares (only 20% 
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photon energy [keV] photon energy [keV] photon energy [keV] 

Fig. 2. looptop (a), footpoint (b) sources and spatially integrated (c) HXR spectra (black) of flare A, 2002 July 23, estimated 
background (gray), and fitting result (red) and its components: thermal (green) and non-thermal (blue). For the looptop spectrum, 
the non-thermal is a thin-target model; for the footpoint and integrated spectra, it is a thick-target model. 




photon energy [keV] photon energy [keV] photon energy [keV] 



Fig. 3. Same as Fig.g] for flare B, 2003 November 02. 



for flare C), as expected. Similar results are found for the full 
(spatially integrated) spectrum. 

The partially ionized or fully neutral chromosphere is 
a more efficient thick- t arget bremsst r ahlung emitt er (e.g. 
lHavakawa & Kitadll95fj: iBrowd 119733: lEmslid fl978l) . which 
can increase HXR flux by a factor of ~ 3 in comparison with 
fully ionized atmosphere. The decrease in the ionization of the 
target reduces the collisional energy loss of the non-thermal elec- 
trons, thus enhancing the efficiency of the HXR bremsstrahlung. 
Applying the thick-target model with a neutral target for the 
footpoint emission, we found the electron rates Nfp /10 35 elec- 
trons i -1 for each flare, respectively, 13 + 2, 60 + 4, 0.4 + 0.1, 
2.0 + 0.5, which are a factor of ~ 2.5 smaller than the values 



found for a fully ionized target. Similar results are found for the 
spatially integrated spectrum. 

Considering the thick-target model with both the albedo 
component and a neutral target, the values for Nfp for each flare 
are 6 + 1, 27 + 2, 0.4 + 0.1, 0.8 + 0.2 in units of 10 35 electrons 
With Nlt found in the previous section, the values for the 
ratio N LT I N FP are 9.3 + 3.0, 10.1 + 2.3, 4.5 + 1 .4, 7.8 +5.2, with 
an average 7.9. 



5. Discussion 

The value Nlt I Nfp > 1 found for the flares suggests either 
(i) the emissivity per unit of length in coronal sources that is 
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Fig. 4. Same as Fig. [2] for flare C, 201 1 February 24. 
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Fig. 5. Same as Fig. |2] for flare D, 201 1 September 24. 



more efficient than free propagation or (ii) additional energy 
losses in the footpoint region. These electron rates also im- 
ply that the main acceleration/injection region is near the loop- 
top, at least for these flares. We would possibly expect to find 
Nlt/Nfp < 1 if the acceleration region is not near the looptop or 
even if the electrons are acce lerated along the loop, due either to 



stochastic electric fields (e. g. Vlahos et al.l20"04T: Turkmani et all 
l2006tlBian~& Brownindl2008l:lGordovskvv et alj|201(fr or to re- 



accel eration by beam-generated Langmuir waves dKontar et aT] 
We now discuss the result of Nlt/Nfp > 1 in terms of 
magnetic trapping, injection properties, and pitch-angle scatter- 
ing. 



5. 1 . Magnetic trapping 

The simplest scenario which can be involved to explain the ex- 
cess of electrons in the loot-top area is magnetic trapping due 
to adiabatic mirroring in a converging magnetic field. The con- 
vergence of the magnetic field, characterized by its mirror ratio 
cr = Bfp/Bit and losscone angle ao = acos(/io) = asin(cr -1 ^ 2 ), 
and the conservation of the magnetic moment of the fast elec- 
trons divide the injected particle population / in a trapped frac- 
tion T, where the pitch-angle a > ao, while the other precipi- 
tating fraction P (with a < ao) escapes the trap and precipitates 
into the dense chromosphere. Defining £ as the fraction of the 
injected population that remains trapped, we have 



/ = T + P, 
T = 



(9) 
(10) 
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Fig. 6. CLEAN maps for the four flares, evidencing thermal emission at 10-11.3 keV (red), footpoints at 61.6-69.5 keV (blue 
contours, at 50% of the maximum of each map), and non-thermal looptop sources (green solid) at 42.8-48.3 keV (50%) for flares A 
and B, 20.7-23.4 keV (40%) for flare C, and 29.8-33.6 keV (50%) for flare D. The countour levels of 30% and 70% at those energies 
(green dashed lines) are also shown. The thick yellow line indicates the length L. Date and integration times are indicated in each 
map. 



P = 

relating this scenario to our results: 

Nfp = P, 
N LT = T + P, 

Nlt _ _J_ 

Nfp W 



(11) where 



(12) 
(13) 

(14) 



N FP 
N lt ' 



(15) 



where the footpoint electron rate Nfp indicates the precipitated 
fraction P, while Nlt indicates the fraction of electrons that pass 
through the looptop, i.e. both the trapped fraction and the elec- 
trons that precipitate directly. 
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Let us consider two typical pitch-angle distributions of en- 
ergetic electrons: (i) isotropic f(fj) = constant and (ii) beamed 
f(p), defined as a Gaussian distribution centred along the loop 
with spread Ap = 0.1. For these two models, one finds the fol- 
lowing relations 



i 



(16) 



(17) 



Eq. [T6land[T7lcan be solved for the losscone value po, and hence 
cr. For the isotropic case, Eq. [T6land[T7lreduce to T/I — po - £, 
and P/I = (1 - po) = (1 - £)■ Solving these equations using the 
observed values of Nu and Nfp, we find the values for £ and cr, 
which are presented in Table [2] Here, we stress that this simple 
model is meant as a way to divide the injected population in 
terms of the losscone, representing an upper limit for the derived 
mirror ratio values, since the accumulation of electrons is not 
considered in the model. 

Under the assumption of an isotropic injection, cr values 
are small (cr > 1 ) and in a greement with values found by 



are small (cr > l ) ana in agreeme nt witn values touna by 
lAschwanden et all (1 19991) and iTomczak & Ciborskil (120071) . In 



the case of the narrow beam, we found that cr > 5, which is 
higher than the range of mirror ratios found by ICosta & Rosall 
(120051) for microwave data, without assuming a particular pitch- 
angle distribution, although their determination is biased by a 
higher weight of the loop-leg source and offers only a lower 
limit of the mirror ratio. When considering our results for the 
case with albedo and neutral target, we found 2.5 < cr < 5 for 
the isotropic case and cr > 17-42, for the beamed injection case. 
As mentioned above, the values for the required mirror ratio are 
just upper bounds. The observations could be explained by lower 
mirror ratios and sustained injection which cause accumulation 
of electrons at the looptop. 

Table 2. Estimated values for trapping fraction £ and mirror ratio 
cr. 



Flare 


A 


B 


C 


D 


ionized target, no albedo 










trapped fraction £ 


0.4 


0.4 


0.3 


0.3 


tr isotropic 


1.2 


1.2 


1.1 


1.1 


tr beam 


6 


6 


5 


5 


neutral target, albedo 










trapped fraction £ 


0.9 


0.9 


0.8 


0.9 


cr isotropic 


5 


5 


2.5 


4 


tr beam 


36 


42 


17 


32 



5.2. Pitch-angle scattering 

The accelerated electrons in the flaring loop are usually assumed 
to be trapped within an adiabatic magnetic trap. Two approaches 
can be used to des cribe their dynamics and spectra, namely 
Coulo mb collisions dMelrose & Browrill 19761: iMacKinnon et alJ 
fl98l and interactions between electrons and waves in the 
plasm a dWentzellll976t iBespalov et aDll987t IStepanov & Tsapl 
120021) . Observed energy-dependent time delays of the 20-200 
keV FTXR emission fav our the Coulomb collision hypothe- 
sis (lAschwandenl 120051) . The time delays of 1 - 10 s for 



smoothed components of the emission follow the law e 3 ^ 2 (e 
is photon energy) and can be explained by Coulomb scat- 
tering of non-thermal el ectrons by dense background plasma 
(lAschwanden et alii 19971) . However, these observed time delays 
can also be explained by a magnetic trap plus pitch- angle scat- 
teri ng by plasma turbulen ce (IStepanov & Tsadl 19991). 

IBespalov et ail d!987l) . followin dKennel & Petschekl i 19661) . 
describe three diffusion regimes of pitch-angle diffusion: weak 
diffusion regime, moderate diffusion regime, and strong diffusion 
regime. The three regimes of pitch-angle diffusion are charac- 
terised by the mean time of pitch-angle scattering tj, flight time 
through the loop t c = L/v, where v is the particle speed, and mir- 
ro r ratio cr = BfplBu/ - weak (t^ > <rt c ), intermediate, or strong 
in iKennel & Petschekl i 1966!) description, (t c < tj < crt c ), and 
strong (f ( / < t c ). In terms of pitch-angle distribution, the weak 
regime generally means that the electrons are weakly (Coulomb) 
scattered in a timescale of the bounce period and that the loss- 
cone is empty. In the intermediate regime, the electrons can 
travel through the loop and are scattered into the losscone, which 
can be considered full, and the pitch-angle distribution is nearly 
isotropic. The strong regime means that a electron can change its 
direction many times during one bounce period. This regime can 
only take place by wave-particle interactions, as a strong colli- 
sional diffusion also means that the fast electrons will quickly 
lose their energy and thermalize. 

To explain the observations in a scenario without magnetic 
trapping, i.e. no coronal convergence of the magnetic field, the 
pitch angle scattering should be strong td/t c ~ {Nlt INfpY 1 , 
with the values of 0.1 < t&jt c < 0.7. Let us consider that the 
HXR-producing electrons are magnetically trapped only at the 
looptop where the emission appears, i.e. the mirror points are 
high in the coronal field. In this case, it is reasonable to expect 
that cr would not be much higher than 2. This scenario can be 
justified because there is no emission detected along the loop 
legs, which would be expected if the electrons were to bounce 
between the mirrorpoints near the footpoints. From our results 
of the trapping analysis in Table [2] we find cr > 2, except 
when an isotropic injection and an ionized target are considered. 
Assuming that at the looptop 1 < cr < 2 but a higher value of cr is 
required to trap the electrons, efficient scattering is needed at the 
looptop to change the pitch-angle distribution and allow the elec- 
trons to be magnetically trapped. Also, the isotropic injection 
implies strong scattering in the acceleration process. Although 
the electron accumulation at the looptop can be explained by a 
trap-plus-precipitation model with modest mirror ratios cr < 2, 
we concluded that the pitch-angle distribution must be wide 
enough to be divided by the losscon e, or else it will r esult i n ei- 
ther a fully trapped (as proposed by iMelnikov et al.l (120021) due 
to an anisotropic pitch-angle distribution orthogonal do the mag- 
netic field, the so-called pancake distribution) or fully precipi- 
tated scenario (as in the classical thick-target model). This re- 
quirement can be fulfilled by moderate/strong pitch-angle scat- 
tering processes. The presence of strong pitch-angle scattering is 
particularly important for stochastic acceleration models, where 
the transport of ele ctrons strongly i nfluences the acceleration 
rate of particles (see lBian et al.ll2012l as a recent review). 

5.3. Instrumental effects 

Because the events we selected for this analysis are GOES M and 
X cla ss flares, it is important to check for pileup (ISmith et al.l 
120021) . This effect, depending on attenuator state, is usually 
stronger in the range 20-50 keV, which is the range where the 
looptop sources were identified. We tested the importance of 
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pileup, comparing the ratio of the corrected and uncorrected 
count rate spectra of each flare. These spectra and their ratio 
are presented in Fig. [7] taken for the time intervals presented in 
TableQ] Presently, it is not possible to correct for pileup in imag- 
ing spectroscopy; however, we can assess its effects in the spa- 
tially integrated spectrum and estimate the errors in the looptop 
and footpoint spectra. Taking the minimum of the ratio between 
corrected and uncorrected count rate spectra as the maximum 
error due to pileup, we find the values of 14% and 20% (both 
around 40 keV) for flares 2002 July 23 and 2003 November 02, 
while errors due to pileup are lower than 2% for flares 2011 
February 24 and 2011 September 24. These errors are much 
smaller than the errors from the fittings propagated to the final 
ratio N LT /N FP : 35%, 24%, 27%, and 64%. Consequently, pileup 
effects do not seem to be significant in our findings. 

6. Summary 

We presented an analysis of the electron distribution in foot- 
points and coronal sources for four well-observed solar flares 
and, for the first time, found the electron rate required to explain 
the coronal and chromospheric X-ray emissions. Using imaging 
spectroscopy techniques, the X-ray spectra of the coronal and 
footpoint sources are obtained and fitted with thin- and thick- 
target models, respectively. The thick-target model directly pro- 
vides the electron rate N F p necessary to explain the observed 
footpoint emission. To obtain the electron rate Nu at the loop- 
top source from the thin-target model, it is also necessary to find 
the length of the source L and the thermal plasma density n (Eq. 
[6]). The length L was measured for each flare in CLEAN maps in 
an energy range where the coronal source was resolved at 50% 
of the maximum of the image and thermal contribution can be 
safely neglected, while the plasma density n was deduced from 
the emission measure EM (from spectral fittings) and source vol- 
ume estimated from ~ 10 keV source size. 

Assuming no albedo and a fully ionized target, we found 
an average NltINfp of 1.6. The ratio NltINfp is further en- 
hanced when co nsidering the Comp ton back-scattered photons, 
i.e. albedo (e.g. iKontar et al.ll2006l) . and a neutral target in the 
chromosphere (e.g. lBrownlll973bl) . Both processes reduce the 
required electron rate to match the observed HXR photon flux, 
effectively reducing Nfp- Under these considerations, the av- 
erage ratio found is ~ 8. The ratio NltINfp> 1 could sug- 
gest that the acceleration/injection region is near the looptop, 
at least for these flares. In previous studies using microwave 
emiss ion maps, homoge neous electron distribution along the 
loop dKundu etalJ 200 1|) or high a ccumulation of electrons at 
the looptop ( Melnikov et al.l 120021) were proposed to explain 
the loop-like morpho logy of the 17 and 34 GHz emission. 
iMelnikov et al.l (120021) concluded that a ratio of electrons at the 
looptop to footpoint in the range ~ 10 - 100 is required to re- 
produce the observations of the four flares considered in their 
study. We note, however, that those results only account for 
the trapped population, without considering the number of pre- 
cipitated electrons required to produce HXR at chromospheric 
footpoints. In fact, their solution requires a pancake pitch-angle 
distribution at the looptop, which would produce very weak 
HXR chromospheric emission (if any at all) due to the lack 
of direct precipitation. Moreover, there is extensive evidence 
of micro wave sources associated with HXR footpoint emis- 
sion (e.g.lKundu etal.ll 19911 iNishio et alll 1 9971: iLee et al.ll2000t 
iKundu et al]l200rj: ISchmahl et alJl2006t iKundu et alJl2009h . in- 
dicating that the trapping conditions vary strongly from flare to 
flare, as one might expect. In fact, evidence for how trapping 



condit ions are different for each flare is presented bv lMeFNikovl 
(1 19941) . showing that the ratio between the number of trapped 
high-energy electrons (i.e. produce microwave emission) and 
precipitated low-energy electrons (i.e. electrons that produce 
HXR emission) can vary by three orders of magnitude, depend- 
ing on the electron lifetime value inside the trap and duration of 
the injection. We note that all those works compared two dif- 
ferent electron energy ranges, while in this paper our results are 
presented for electrons in the deka-keV energy range. 

We considered a collisionless magnetic trapping to explain 
the electron rates found (for fully ionized target and no albedo) 
and found that the required mirror ratios are in the range 1 . 1 < 
cr < 6, while for a neutral chromospheric target and account- 
ing for albedo photons, 2.5 < cr < 17 - 42. These values reflect 
an upper boundary for the mirror ratio, as our simple model do es 
not account for electron accumulation over a sustained injection. 
Therefore, the observations can possibly be explained consider- 
ing a trap-plus-precipitation scenario with modest mirror ratio 
values. However, we note that the pitch-angle distribution must 
be wide enough to be divided into trapped and precipitated frac- 
tions by the trap losscone. This requirement indicates that mod- 
erate/strong pitch-angle scattering must take place in the accel- 
eration site or during the transport along the loop. 

If we consider a scenario without significative convergence 
of the coronal field towards the chromosphere, i.e. no efficient 
magnetic trapping, the pitch-angle scattering time required to 
explain the observations appears to be in the range 0.1 - 0.7 of 
the time required for an electron to cross the length of the loop. 
The observations do not support very strong scattering with val- 
ues smaller than 0.1, which is sometimes required for efficient 
stochastic acceleration. This conclusion indicates that the ob- 
served coronal source is not the acceleration region, which could 
be in a smaller or larger portion of the loop. Hence the coro- 
nal HXR source is mainly due to transport processes. Equally, 
it can also mean that the stochastic acceleration models with 
strong pitch-angle scattering are not the dominant mechanism 
for electron acceleration at these energies. Conclusions in favour 
of the tr ap-plus-precipitation model are given by other stud - 
ies (e.g. lAschwanden et al.l 119991: iTomczak & Ciborskill2007l) . 
Nevertheless, the effect of moderate/strong scattering can possi- 
bly explain the formation of the HXR looptop source or at least 
enhance the effect of the magnetic trapping. It is our impression 
that the scattering cannot be neglected and should be investigated 
with the trap-plus-precipitation model for the deka-keV energy 
range in more detail. 

In conclusion, our results suggest that the accelerated elec- 
trons must be subject to magnetic trapping and/or moderate or 
strong pitch-angle scattering, keeping a fraction of the popu- 
lation trapped inside the coronal loops. This is in full agree- 
ment with the trap-plu s-precipitation models first proposed by 
iTakakura & Kail d 19661) and further developed by many other au- 
thors. 
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